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ABSTRACT 


Prequency domain methods are employed for the design and 
analysis of anti-roll tank systems. A waveslope spectrum of 
the ocean is developed based upon the trochoidal model of 
ocean waves and the Pierson-Moskowiltz ocean wave power 
Spectra. The response of a typical ship to the waveslope 
spectrum is studied with the aid of a numerical example. 

FORTRAN subroutines for the design of passive and 
pump-activated U-tube anti-roll tanks are developed, and the 
effects of a tank installation on the ship parameters is 
observed. Many controller design methods are investigated, 
and final design is accomplished by exhaustive search of 
closed-loop frequency responses. Frequency responses of 
passive and activated tank stabilized ships are compared 
With that of an unmodified Ship to Lllustrate the 
effectiveness of design. The effects of vessel displacement 
Variations due to consumption of fuel and stores on tank 


effectiveness and controller design are discussed. 
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I. FORMULATION OF THE PROBLEM 


A ship on the ocean 1S in constant motion. Ideally, the 
only motions would be translation in the forward direction 
and rotation about the vertical axis, associated with 
propelling and maneuvering the ship. Realistically, there 
are perturbation motions in addition to the desired motions 
in the form of translations along and rotations about the 
longitudinal, transverse, and vertical axes of the ship. 
Th2 purpos2= of this thesis was td study the rotation about 
the longitudinal axis, Known as roll, and arrive at a scheme 
poe COntrolling sige A procedure for designing the 
installation of a U-tube anti-roll tank system was 
consolidated from several sources | Chadwick, 1950 | [ Chadwick 
wee Xlotter, 1953 || Blagoveshchensky, 1962 || Vachanaratana, 
1973 || Ziegler, 1975] and presented in the form of a FORTRAN 
Subroutine. Many controller design methods were 
investigated, and the final design was accomplished by 


exhaustive search of closed-loop frequency responses. 


Bee WAVESLIPE SPECTRUM OF THE OCEAN. 


Before a meaningful analysis of ship motions in water 
could be conducted, the motions of the water itself had to 
be understood. Because a Ship tends to ride normal to the 
surface of the water, the most Significant physical property 
of waves with respect to ship motions is the waveslope, y. 
At the crest and in the trough of a wave, the waveslope is 
zero, while at some point nearly midway between those 
extremes, the waveslope is maxinun. 

Several theories of wave motion have been advanced, 


including the simple harmonic wave; the Stoke's Wave, which 





is a harnonic wave with a second harmonic added; and the 
trochoidal wave. Geometrically, a trochoidal curve is 
constructed by rolling a circle of radius R under a 
horizontal line, and tracing the locus of a point a radial 
distance c<R from the center of the circle. The trochoidal 
curve may be steep or shallow depending on the ratio of r:R, 
and corresponds very closely to actual observed ocean waves, 
where 27 R = L = wavelength, and 2r= h = trough-to-crest 
wave height. 

Ocean waves are primarily caused by the transfer of 
kinetic exergy from the wind, which is the result of the 
transfer 9f thermal energy in weather phenomena. As_ the 
wind blows over a flat surface of water, the water particles 
at the air-sea interface are excited by frictional forces, 
and tend to move downwind. The freedom of motion of these 
particles is limited by the fluidic properties of water, and 
thus they tend to accumulate in successive ridges, or 
ripples, on the sea surface. Once these ridges are formed, 
further wave buildup continues by a different mechanisn. 
On the upwind side of a ridge, a small local high pressure 
region develops which exerts a downward force on the water 
surface, forming the wave trough. On the downwind side of a 
ridge, a corresponding small local low pressure region 
develops which draws the water surface upward, forming the 
wave crest. The wave continues to grow as long as the 
wind continues to supply energy in excess of that consumed 
by the internal friction of the water in the wave, which it 
does until the wave velocity equals the wind velocity. 

Three key factors govern the development of wave 
systems. They are the initial difference between the 
velocities of the wind and any existing wave systems; the 
length of time (duration) that the wind continues at a given 
Mewocity; and the fetch, or distance of relatively uniforo 
depth water over which the wave will run before being 





altered by a change of depth or destroyed by encountering a 
land mass. (The passage of a vessel through an ocean wave 
system does little to alter the growth or shape of the 
system.) 

The rate of increasing wave height is a function of the 
initial difference between the velocities of the wind and 
the existing wave system. A sudden change of wind velocity, 
such as that associated with the passage of a front, can 
cause the waves to increase at a rate of one or two feet per 
Minute. There is a physical limit to the increase in wave 
height, aad it is defined in terms of the wave height to 
wave length ratio, h:L=1:7. Near this ratio limit, excess 
Wind eénecgy blows off the tops of the waves, forming 
turbulence known as "white caps." 

The rate of increasing wave length generally lags behind 
the rate of increasing wave height. Thus in regions of 
short fetch, a sudden change of wind velocity may never 
produce waves with long wavelengths due to the limiting 
ratio. In contrast, on the open ocean where there is a long 
fetch, a change of wind velocity can produce long waves of 
relatively low height. For waves of length L>100 feet, there 
is a diminishing likelihood of ever approaching the maxinun 
value of the ratio. [In practice, for wavelengths between 
500 and 300 feet, the ratio rarely exceeds 1:20; for the 
longer ocean waves (of length over 1000 feet), the ratio 
drops to 1:50. 


The above discussion of wave development only hints at 
the complexity of attempting to define the waveslope 
spectrum of the ocean. Conferences have been held, and many 
books have been written on the subject of ocean spectra. 
| conference, Eason, Haeyland, 1961 | | Inoue, 1966 | 
| Comstock, 1967 | | Gillaer, 1969 | | Ziegler, 1975 | Many ocean 
Spectra plot relative wave energy as a function of frequency 
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with "sea state" as a parameter. Others define wind 
velocity and duration as parameters. Stiil others define 
"Significant wave height" as a parameter, with many 
definitions of "significant wave height" itself. Spectra 
have been constructed for given locations citing mean wave 
height, or wean squared wave height, or RMS wave height as 
pactameters. 

Two other factors which affect the effective ocean 
Spectra are vessel speed and angle of encounter. To 
illustrate these effects on the ocean spectrum, consider a 
Ship proceeding with velocity identical to that of a 
Single-frequency wav2 train. If the ship heads directly 
into the approaching wave train, the apparent frequency of 
the wave doubles, but tne roll coupling with the ship goes 
to zero. When the ship travels parallel to the wave train, 
which 1s to say the ship 1s ina beam sea, the apparent 
frequency of the wave is the same as the actual frequency of 
the wave, and the roll coupling is maxinun. When the ship 
heads in the same direction as the wave train ("following 
seas" coniition), the apparent frequency of the wave drops 
to zero, as does the roll coupling. If B denotes the ship's 
angie of encounter with the wave train, where B=0 fora 
following sea, B=90° for a beam sea, and B=180 for a head 
sea, then the apparent waveslope is related «to the actual 
waveslope, and the apparent frequency of the sea is related 
to the actual frequency by the expressions | Ziegler, 1975 |: 


¥(B) = ¥ sin B (1) 
maxX 

w(B) = w (1.- w uu. /g) cos B (2) 
sea sea ship 

Wes cOnwanrda velocity of the ship 

ship 

g = acceleration due to gravity 


Since the ultimate reason for this study was to aid in 


nies 





the desigi of a general-purpose anti-roll stabilizer, the 
ocean spectrum chosen was a wave height power density 
Spectrum, with sea state and significant wave height (1/3 
highest waves) aS parameters. The waveslope spectrum was 
then computed on the basis of the trochoidal model of the 


ocean, where 


r = 1/£ = 27r/w_ seconds (3) 
2 

L = 5.118 T feet (4) 

v = (h/L) x 57.3 degrees (first SRE DORT NETO 

max (5 


Be. SHIP ROLLING IN WAVES. 

A ship on the ocean moves, in general, with six degrees 
of freedom relative to the coordinate axis system. By 
Convention, the positive x-direction 1s forward along the 
longitudinal centerline of @ehe Ship the positive 
y-direction is horizontally outward on the starboard bean, 
and the positive z-direction is vertically downward. The 
center of the coordinate axis system is located at the 
vessel center of gravity. Rotation about the x-axis is 
Known as roll, with positive roll angle, ¢, defined as an 
inclination to starboard. Roll can be considered separately 
from the other degrees of freedom, which reduces the study 
of ship motions to a two-dimensional problem with one degree 


of freedon. 
1. Definitions and Concepts of Naval Architecture. 


The following definitions and concepts of naval 
architecture | Comstock, 1967 | |Gillmer, 1969 | are given for 
the benefit of the reader, and are illustrated in Figure 
1-1. 

a. Displacement. 

The displacement of a vessel, A, is defined as 
the weight of the volume of water displaced by the vessel. 


Ve 





WL WL 





eee— Center of Gravity 

meee Center of Buoyancy 

feee— Metacenter 

GM = Metacenter Height 

GZ = Righting Arm 

@ = Ship Roll Angle Relative to Earth Reference 
WL = Water Line 


Figure l-l. Definitions and Concepts of Naval Architecture 
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Due to the Principle of Archimedes, the displacement equals 
the weight of the vessel itself. Displacement is always 
given in long tons, where one long ton equals 2240 U.S. 
pounds. The weight of the water displaced equals the weight 
of the vessel, independent of the density of the water. The 
volume of water displaced, however, is a function of the 
density of the water. For purposes of most calculations, 
the density of fresh water is given as 36 cubic feet per 
ton, whereas the density of salt water is given as 35 cubic 
feet per ton. 
Dee Genter of Gravity. 

The center of gravity, G, is defined as the 
point at which the summation of the moments of all the 
weights ina vessel with reference to any axis passing 
through that point equals zero. The center of gravity of 
most warships of usual form is generally located near the 
designer's waterline, in or near the plane of the midship 
section. For purposes of buoyancy and stability analysis, 
the entire weight of a vessel is considered to be 
concentrated at, and the force of the weight acts vertically 
downward through, the center of gravity. 

c. Center of Buoyancy. 

The center of buoyancy, 8B, is defined as the 
geometric center of the underwater form of the vessel. The 
force of buoyancy, equal in magnitude to the force of the 
displacement of the vessel, acts vertically upward through 
the center of buoyancy. When a vessel is in static 
equilibrium, the center of buoyancy is located directly 
under the center of gravity, such that the force vectors of 
displacement and buoyancy are equal and opposite. Due to 
this definition of the center of buoyancy, as the condition 
of a vessel varies from that of static equilibrium, the 
location of B varies also. 
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d. Metacenter. 

As a ship inclines, the shape of the underwater 
form chanjy2s, and the center of buoyancy migrates to remain 
at the geometric center. The migrations form a locus which 
is dependent upon and directly related to the shape of the 
Ship's hull. This locus is generally elliptic and as such 
Cannot have a fixed center, but for an infinitesinal 
increment of the locus the instantaneous center, calleda 
metacentec, can be located. The locus of metacenters is 
geometrically an evolute, andis called the metacentric. 
For small angles of inclination, up to 10° for most ships, 
the metacenters are co-located at a point known as the 
Metacentec, 4M. 

e. Metacentric Height. 

When a ship is upright, that is, when the angle 
of inclination = 0 , the centers M, G, and Blie ina 
straight vertical line over the keel. The distance between 
Gand Mis known as the metacentric height, GH. If Mis 
above G, 3M is defined as being positive. 

£f. Righting Aram. 

When a ship is inclined, the center of buoyancy 
1s no longer on the original vertical centerline. A point Z 
is described by the intersection of a line drawn vertically 
upward from Bwith a line drawn horizontally away from Gc. 
The distance GZ is known as the righting arm, and for small 
angles of inclination 9, 

GZ = GM sing = GM > (6) 
g. Righting moment. 

The righting moment, M,, is the product of the 
force (displacement) and the righting arn: 

MM. = AGZ= AGM (7) 
h. Metacentric Height and Stability. 

The metacentric height is a measure of the 
Statical stability of a vessel. A large positive GM 
corresponds with a pair of high-frequency complex conjugate 
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poles in the left half of the s-plane, far removed from the 
jw -axis. A small positive GM corresponds with a pair of 
low-frequency complex conjugate poles in the left half of 
the s-plane, close to the j -axis. If GM = OQ, then there 
are two poles at the origin of the s-plane. A negative GM 
puts poles in the right half of the s-plane. 


Zee Equations of Motion. 


as Natural Period of Roll. 

In still water, a ship can be made to roll by 
applying and removing external heeling moments. In the 
inclined position, the righting noment is equal and oppcesite 
to the heeling moment, and the potential energy of the ship 
is equal to the work done by the external force causing the 
inclination. When the external moment is removed, the 
righting moment produces rotation of the ship towards the 
mepagot o0Sition. The potential enerqy is converted to 
Kinetic energy such that when the ship is in the upright 
position, assuming unresisted roll, all the potential energy 
has been converted to kinetic energy. As a result, the ship 
Will continue to rotate past the upright position untii all 
kinetic energy has been reconverted to potential energy, at 
the opposite limit of roll. Assuming no loss due to 
Mmenection, the ship would continue to oscillate, or roll, 
indefinitely with constant amplitude. The period of these 
oscillations can be found from simple harmonic motion 
relationships to be | Comstock, 1967 | 


T = (27k) /A/gGM 
7 (1.108k) 7. /Gu 
=~ (0.44b)//GM- (8) 
where 
g = acceleration due to gravity, ft/sec 
k = radius of gyration of mass of ship, feet 
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b = extreme beam of ship, feet 


b. ROl1l by Wave Action. 

The cause of ship rolling in a seaway is 
primarily the unbalanced moments which result from the 
center of buoyancy shifting in response to the change of the 
waterplane, and thus the shape of the underwater portion of 
the hull, as waves pass under the ship. A ship will usually 
attempt to ride normal to the surface of a wave. A ship 
also moves in a circular orbit as do the particles of water 
in a wave. The orbits induce centrifugal forces, which are 
opposed by the dynamic forces in water. 

Rolling in waves is a forced oscillation, and the ship 
will tend to oscillate at or near the frequency of the 
waves. If the waves are not uniform, which is usually the 
case, aship will attempt to oscillate with its natural 
frequency. The overall roll, therefore, is a combination of 
the ship's natural and the wave's periods, with the wave 
influence being generally more predominant. 

Cc. Equations of Hotion. 

A ship can be modeled as a second-order, 
Spring-mass-damper systen. | Chadwick, 1950 ] | Comstock, 1967 | 
The spring constant, or static righting coefficient, was 
discussed earlier in this paper in the definition of the 
righting moment. The static righting coefficient is 

K =AGM (9) 
S 

The inertia term, or roll inertia coefficient, is the 
product of the mass of the ship and the radius of gyration 
of the mass of the ship squared: 


2 
ue = (A/g) k (10) 


The radius of gyration is found by timing the natual 
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period of roll in protected (still) water, and massaging the 
equation for roll period as follows: 


kK = (Tg (GGH )/2m 


= (Ty/GM )/1. 108 (11) 
The coll damping coefficient is derived from the 


observation that the damping ratio of a typical ship is 
approximately 
t = 0.1/7 = 0.0318 (12) 
ship 
By working backwards through the second-order systen 
relationships, the roll damping coefficient was found to be 


2 
B = 0.1 (1,4GM) /fr (13) 


Ss > 

On the forcing function side of the equation of motion, 

following the classical usage of Froude and Krylov, 

| Chadwick and Klotter, ee the ship-sea coupling 

coefficient was considered to be approximately equal to the 

Seatic cighting coefficient. Therefore, the equation of 
motion for a ship in a seaway is: 


JI¢ +BhO +K GH =FK YW (14) 
S S S 

Dividing through by the roll inertia coefficient and 
performing the Laplace Transform operation, the system 


transfer function was obtained: 


2 
[s + (B_/d_)s + (K_/J_)] @(s) = (K_/J_) w(s) (15) 
Ss Ss S Ss S A) 
2 2 2 
[s + tw )s + # |] o(s) = © w(s) (16) 
s Ss Ss ) 
2 
&(s)_ w 
2 
Y (Ss) Ss + (2 w)s + w 
s Ss Ss 


ES 





K /J 
re 2) (17) 


Z 
s + (B /J )s * (K /J ) 
s s s <« 


It was useful in the analysis to normalize the equations 
such that the ship's natural frequency was unity. This was 
accomplished in the following manner: 

Starting with the general equation of motion, divide 


through by the static righting coefficient: 


(J /K )od+ (B SK )db+t Go zY (18) 
S Ss = S 
meme = oy £ = (/K /d jt 
S Ss S 
Define ado /dt = ¢ 
dod /dt' = ¢' 
2 2 
dgod/dtt = >" 


Then $= (/K /J ) 
s 5s 
= (K /J _) >" 
Ss s 
The eyuation was then rewritten as 
o"+ (B/S KI ) e+ oO = ¥ (19) 
AS) Sos 


The LaPlace Transform operation was performed, and the 


System transfer function (normalized) was obtained: 


2 
E + (B/S KJ )st¢ 1] ® (s) = y(s) (20) 
S S $s 
r ees ae 
VY (Ss) s + (B/KJI)s +1 
S Ss Ss 


3. Fcequency Response. 


In order to study the frequency response of a vessel 
to the waveslope spectrum of the ocean, it was necessary to 


introduce a numerical example. The example chosen was 4 
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well used example | BLlagoveshchensky, 1962], which was: 


L = 49.1 meters (161.0 feet) 

b = $.0 meters ( 29.5 feet) 

&O = 936.0 metric tons (921.2 long tons) 
GM = 0.73 neter ( 2.4 feet) 

Ty = 8.63 seconds 


The radius of gyration of the mass of the ship, and the 


coefficients of the eguation of motion were computed as 


follows: 
k = (8.63.\/(9. 81) (0.73) )/2r = 3.68 meters 
Z. 2 
J = (936/9.81) (3.68) = 1289.021 tons-meters-sec 
S 


Z 
B = (0.1) (8.63) (936) (0.73) /7 =59.7 tons-meters-sec 
Ss 


K = (936) (0.73) = 683.28 tons-meters 
Ss 


These values were then substituted into the transfer 
Function, the roots of the characteristic equation were 


found, and the frequency response was determined. 


@(s) = Ons 3 


eee 22 
wy (S) s + 0.046s + 0.53 ae 


roots (poles) at See eUZ ee ) 0.73 
natural frequency w = 0.73 
S 


damping ratio cr = 0.032 
S 


The frequency response is graphically illustrated in 
Pigure 1-2. The Digital Simulation Language (DSL/360) 
program which generated tne plot is listed in Appendix A. 
The zero decibel low frequency response demonstrates the 
natural tendency of a ship to ride normal to the surface of 
a wave. The resonant peak which occurs at 0.724 rad/sec has 
Megan Of 23.9 db, or a Wagnification factor of 15.6. 


According to the waveslope spectrum of the ocean, tabulated 
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in Appendix E, for a sea of 15 foot significant wave height 
(s2?a state 6) the maximum waveslope angle at 0.724 rad/sec 
would be 2.23 degrees. Therefore, the ship would experience 
@ maximum roll angle of 34.8 degrees, relative to the 
horizon. This severity of roll causes crew discomfort, and 
limits the ability of the ship to conduct operations, such 
as gunfire support, flight Operations, or small boat 


operations. 
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Piste LADEeGEGALLON BY PASSIVE ANTI-ROLL TANKS 


A. DEVELOPMENT OF THE EQUATIONS OF MOTION 


In order to apply servo theory to the problem of ship 
stabilization, the physical system was. assumed to be linear, 
with lumped parameters and constant coefficients. It was 
also assum2d that the rolling degree of freedom could be 
considered independently of the other degrees of freedon, 
such that the ship oscillates as if about a fixed center of 
rotation. The problem thus became a two-body problem, the 
Ship and the tank-fluid, each with one degree of freedon. 
The tank-fluid was idealized as a one-dimensional filament, 
With provision for varying cross-section but with constant 
velocity across each cross~section. These assumptions 
reguired small angles of ship roll but allowed for moderate 
angles of tank-fluid level, both being reasonable in a 
ao lizel ship. 


1. Double Pendulun Analogy. 


As discussed earlier in this paper, a ship at sea 
behaves like a damped oscillator, and was modeled as a 
second-order spring-mass-damper system. Another appropriate 
moiel which was considered was the pendulum model [ Chadwick, 
1950 | [ Chadwick and Kotter, 1953 |, SaecnecdesDring constant, 
an inertia constant, and a damping constant. 

The stabilizer chosen was a U-tube anti-rolling tank 
System. Th2 water or other fluid in the U-tube acts asa 
nydraulic pendulum. Since the tank was to be contained 
Within and wholly supported by the ship, the ship-tank 
System that resulted was a pendulum supporting a pendulun, 
Or a double pendulum. In general, this kind of syst2m leads 
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to two second-order differential equations with coupling 


terms. 
2. Nomenclature. 


There are two types of nomenclature needed here. 
These are the parameters and variables of the differential 
equations, andthe dummy variables and parameters which 
characterize the tank geometry. 

ae Differential Equations. 

The differential equations are second-order in 
nature, and have forcing functions. The variables and 
parameters are as follows: 

Pepeident miactaples of the forcing Side: 
Y = instantaneous effective waveslope 
Paramaters of the forcing variables: 


K = ship-sea coupling coefficient 


Ss 


Torques produced on the forcing side: 


K wp = disturbance torgue due to waves 
ss 


Dependent variables of the homogeneous side: 


Q 
0 


Ship roll angle relative to the horizon 


tank-fluid level angle relative to the ship 
Parameters of the homogeneous variables: 


Ship self-terms; 


J' = roll inertia coef., unmodified ship (no tanks) 
B' = damping coef. of unmodified ship 

S 

os Stdtie Eeghtimg coef, Ofeuntodified ship 

S 
J = roll inertia coef., modified ship (with tanks) 
e = damping coef. of modified ship 

K = static righting coef. of modified ship 


Vj 
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Tank self-tearms: 


7 = equivalent inertia coef. of tank-fluid 

7 = skin-friction damping coef. of tank~fluid 
a = overall damping coef. of tank-fluid 

- = equivalent static coef. of tank-fluid 


Matual coupling terms: 
J ee mutual-coupling inertia coefficient 


B BS mutual-coupling damping coefficient 


K os Mutual-coupling static coefficient 


It is a feature of the double pendulum that 


K = K 
st ite 
B = Q 
st 
As discussed earlier in this paper 
K = K 
3s S 


db. Tank Geometry. 
The notation used in the derivation which 
depends upon tank geometry is shown in Figure 2-1, and is 


given below. 


S a generalized coordinate, the curvilinear 

distance along the fluid trajectory. s=0 at 

the tank-fluid surface in the starboard tank, 

and s=S at the tank-fluid surface in the port 

tank, when @ =0. 

R = lever arm of tank, the horizontal distance 
from the vertical centerline of either side 


tank to the vertical centerline of the Ship. 


r(s) the radial distance from the virtual center 
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Figure 2.1. The Dependent Variables and Nomenclature of Tank 
Geometry 
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of rotation to any Poun-es on the fluid 


ehagectrOly. 
p = density of fluid in tanks. 
D(s) = the perpendicular distance from the virtual 


center of rotation to tangent at any point on 


trajectory. 
a(s) = angle between r(s) and D(s). 
y¥(s) = angle between r(s) and deck of ship. 
A = tank cross-section at s=0 and s=S. 
a(s) = tank cross-section at any point s. 
g = acceleration of gravity. 


3. Euler-LaGrange Equations of Motion. 


Ti2 two degree of freedom system of ship and 
tank-fluid was best treated by means of the Euler-LaGrange 
equations in generalized coordinates. The generalized 
equation is | Chadwick and Klotter, 1953 | 


a fat at au ap 
tog taaget oq + 9q =Q . (23) 
ae 1 1 ale le 





where Q; is the so-called generalized non-conservative force 
in the ith equation, and i=1,2. The generalized force is 
equal to the rate of doing virtual work against the external 
mon-conservative forces when q; is given a virtual 


displacement. The parameters of the Euler-LaGrange equation 


are: 
- ; 
=> ny © 72 = kinetic energy (24) 
=») J tax = eae eee energy 
(due to gravity) (25) 
2 
D = By /2 = dissipative ener 26 
2 4 (due to skin ame oon etc.) i 
£ = conservative force 
B = damping parameter 


ae 





4. Calculation of the Energy Functions. 


For each type of energy (kinetic, conservative 
potential, and dissipative), the contribution of the ship 
alone andthe contribution of the tank-fluid alone was 
calculated, then the contributions were added to arrive at 
the total energy. In the regions where the fluid surfaces 
travel, namely the vertical legs of the U-tube, it was 
assumed that the tanks were straight, parallel, and of 
uniform cross-section. This meant that a(s) and D(s) were 
constants in these regions, thus for any expressions which 
depended only on the variables a(s) and D(s), the limits of 
integration were fixed at zero and S, without approximation. 


Calculation of kinetic energy: 


S+R tané 5 5 
T = foa(sy)/6R aAl+| dr(s)| + 2/6R alld d(s)| (ds 
fluid 2 a (s) a (s) 
R tané (27) 
-2 
A. = J' d /2 (28) 
ship S 
a = i at (29) 
total f bond ship 


Calculation of potential anergy: 


S 
U = |p A rf 2 S) r(s) siny dsii1 - cos @¢ 
fluid : a || | 
5 (30) 
+ pgAR 2 tan@ sind@- tan @ cos? 
2 
ut =  K'd /2 (31) 
ship S 
; pee EU. (32) 
2 otal a1 d ship 
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Calculation of dissipative energy: 


D = B'@ /2 33 

fluid t ‘ a 
°2 

De = B'd /2 (34) 

ship S 

D = D Loot ; (35) 

total Elba d ship 


eens Equations of Motion. 


The equations of motion of the ship and tank systen 


cone from equation (23), =with ds = @?, and q = 6. 
Substituting in equations (27) through (35) and simplifying, 
and letting Q, = Kv and Q, = 0, the equations of motion are 
Ship: (J $+B 6 +K 6) + (J O6+K 6) = KW (36) 

S Ss S st ie Ss 
Tank: (J $+K 6) + (J 6 +B 6 +K 8) = 0 (37) 

st ie te e c 

Obtaining the System transfer Mine 2on {Lov the 


equations of motion was somewhat more complicated than it 
was for the ship without tanks. First, the equations were 
divided through by the respective inertia coeficients, and 


the LaPlace Transform operation was performed: 


Ship: [s+ (3 aK 73 yJecsy+ [cg /J)s+(K /3 ) ]e(s) 
s s s s st s t s 
= (K /J ) w(s) (38) 
Ss Ss 
Tank: Ke / J eae /3 ) |e (sy +[s°+ (B /J jpst(K /J ) |ocs) 
st ¢ tt tt tt 
= 0 (39) 


@(sS) was isolated by applying Cramer's Rule, and the 


transfer function was obtained: 
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J J -J 
Ss Ss 


st ts < 
B K #K B K K -K 
s Ss S 
s + 2 (40) 
el J J ~-J 
s t st s t st 
b ) (K /J ) “+B /3J )st+(K /J ) 
S s S 
_ aS = s s ee OE tat (41) 
Y (5) A 


6. Definitions of Parameters. 


The parameters of the equations of motion are given 
in terms of previously defined parameters and dummy 


Variables as follows: 


Ship self-terms: 


S 
2D 2 
J = J' + paAR 
: ee {aoa ds (42) 
0 
B = B (43) 
S Ss 
S 
Kk = Ki + pgaR fp [rz (s) 7% | siny ds (44) 
0 
Tank self-ternms: 
’ S 
i, = paAR f A/a(s) ds (45) 
0 
B = B? 46 
t = a 
2 
Ko = 2Pg9AR (47) 
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Mutual-coupling terms: 


S 
2 

y es p AR f D(s)/R ds (48) 
~ 0 
B = Q 4g 
st oe) 
K = K (50) 
st t 


Peo NALYSIS OF THE EQUATIONS OF MOTICN. 


Having arrived at a set of 2quations of motion, it was 
necessary to conduct a brief analysis of the equations to 
give them physical significance. 


1. The Self Terms. 


Looking first at the self-terms of the homogeneous 
equations, it was observed that the ship acts as a simple 
pendulum if the tank-fluid does not move, and that the 
tank-fluid acts as a simple pendulum if the ship is 
prevented from moving. Each of these pendulums has its own 
matural frequency. These two frequencies are important 


factors in the stabilization problen. 
2. The Mutual Ternaus. 


There is a third critical frequency associated with 
the ship-tank system, called the secondary resonance 
freguency, which is due to the mutual-coupling terms. The 
mathematical significance of the effect became apparent when 
the ship equation, equation (36), was rewritten with the 
Mutual tecms on the right side as follows: 


JI O+*B OK P= KY -(T FF +K 8B) (51) 
Ss Ss Ss Ss st uv 


out 





This 2quation shows the ship as a simple pendulum, acted 
upon by a net torque, which is the torque due to waves minus 
the torque due to tank-fluid motion. If the tank-fluid were 
made to oscillate at the secondary resonance frequency, the 
tank-fluid torque would be equal to the wave torque, and 


thus the ship would be influenced by zero net torque. 


C. ANTI-ROLL TANK DESIGN CRITERIA. 


The easiest tank geometry to design and build is a right 
rectangular U-tube, with the vertical legs having constant 
cross-sectional area, and the horizontal cross-leg having 
constant but smaller cross-sectional area. Given this tank 


geometry, 
a(vertical leg) = A (52) 
a(horizontal leg) = a ; 5 
= (8 ee - h a (53) 
D(vertical leg) = (54) 
D(horizontal 12g) = (55) 


The length of the trajectory of the fluid in th2 tank, 

S, is dependent upon the lever arm of the tank, R, and the 

height of the fluid column in a vertical leg of the U-tube 
when @ = Oo - les 

S = 2(R + h) (56) 

For maximum effect, the tank system must span _ the 

extreme beam of the Ship. Therefore the lever arm of the 

tank is half of the beam, b, diminished by the width of a 
vertical leg, w, 

R = (b-w) /2 (57) 

The relationships derived fron the equations of motion 


which define the tank dimensions and parameters are: 


er 





K f/J = A 58 
: (58) 


Ss 
~ We 

K, 75. = 
f A/a(s) ds (59) 
0 

J /35 = X 

st ite 
X = 


S 
fo a)/(R A) ds 
0 


Zee e 2 (a7 A) (60) 


meme rist Latio, \, E@tlects the free surface effect of 
the tank-fluid, and should be a small number. An acceptable 
value was found to be Q.05 | Vachanaratana, 1973 |. 

To be effective as a stabilizer, the anti-roll tank 
system must be "tuned" to the natural frequency of the ship 


in which it is installed. Therefore, 


2 2 
GW = WwW 
te 
K /J = K JSJ (61) 
te ie s s 


The third ratio, X, is dependent upon the placement of 
the tanks in the ship relative to the center of rotation. 
This ratio can be forced to zero (J, forced to zero) by 
placing the horizontal leg of the U-tube above the center of 
rotation by the amount 


-D =h A/a (6 2) 


There are no hard and fast rules for determining a value 
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of h, except that the tanks must fit within the hull of the 
ship, and the fluid in the tanks should not be permitted to 
either contact the top or completely drain out of the bottom 
of a vertical leg of the tank at the maximum expected 
tank-fluid level angle, Onay 

Since J, depends upon tank geometry, the approximation 


must be made that Ue = Jy , for purposes of initiating 
design. Then 
K = i J! 
e Ss 
2 
= 2pgGaAR 
2 
= pglw (b-w) /2 (63) 


Aa reasonable value of tank width should be chosen. Then 
Z 
1 = Z SS f pq ww (b-w) (64) 


Making the approximation that the natural frequency of 
the modified ship is the same as the natural frequency of 
tha ship without tanks, and picking a value for h, equation 
(53) is solved for a, and then equation (62) is tested for 
feasibility. If the value of D so computed cannot be 
physically realized, then a realizable value of D_ must be 
chosen. This merely means that Jey Will not be identically 
zero, and the overall analysis will take more care. 

The tank damping parameter is derived by working 
backwards through the second-order relationships for the 
hydraulic pendulum, noting that the damping ratio for the 
tank-fluid is approximately equal to the damping ratio of 
the (unmodified) ship: 


cn = 0.1/7 (65) 


B = 0.2 /K J / 66 
: V i, © (66) 


fe 
Once the tank system parameters have been determined, 
the modified ship parameters are found. Over the range s=0 
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to s=h, r! m su (D+h-s) e and xr sinY= Dth-s; over the 


range s=h to s=htR, rf = (R+h-s)° + me peand © Sin Y = 0. 
Z 
J = J' + 2pAR (I + I) (67) 
S S 1 2 
Z 2 2 3 2 
r = [ (3R +2D )h + 2Dh +h |z [ 32 ] (68) 
2 2 
= = a (R +3D )/( 3AR ) (69) 
K =  Cteec pg A Rk (I + I ) (70) 
S$ S 3 4 
2 
* = -(2D + h )/( 2 ) (71) 
=n = Da/Aa (72) 


The design steps should then be retraced, removing the 
approximations. Normally, after two or three iterations of 
the design procedure, the calculated values converge on the 
actual values, anda final design has been achieved. The 
FORTRAN subroutine used in this iterative design procedure 


is included in Appendix B to this thesis. 


Bee CE OUENCY RESPONSE OF PASSIVE SHIP-TANK SYSTEM. 


The same numerical example used in the frequency 
response study of a vessel on th? ocean was used for the 
design example and frequency response study of a ship with 
passive anti-roll tanks to the waveslope spectrum of the 
ocean: 


Jt = 1289.021 metric tons-meters-sec 
3 

B! = 59.746 metric tons-meters-sec 
$ 

K' = 683.280 metric tons-meters 


The following arbitrary choices were made: 
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= 1.0 
= 1.3 
= 0.05 


oT 3 FF & 
| 


3 = 1.025 


Then 


A 
tt 


meter 


meter 


3 
m2tric tons/meter (salt water used) 


(0.05) (1289. 021) 


= 64.451 metric tons-meter 


Z 
Ll = (2) (64.451) / (1.025) (1) (9-1) 


= 1.965 


A =) Weel 


= 1.965 


meter 


meter 


a = Che ee A (63 3 oes e tT 


2 
0.457 meter 


Testing equation (62) for feasibility, 


=D 
eo 


Since this tank 


in the ship being 


above the center of 


(1.3) (1.965) / 0.457 


meter 
placement would probably not be possible 
studied, a more realistic value, still 


rotation, was arbitrarily chosen: 


meter 


Then X is non-zero, as indicated. 


J /3J = X 


st t 


2 [1.3 - (0.457) (1.965) |] 


= 2.135 
To summarize, the tank dimensions found thus far were 
w = 1.0 meter 
Ll = 1.965 meter 
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2 
1.965 meter 


A = 
aoe = O. 457 Reree- 
R = 4.0 meter 
D = -1.0 meter 
Given these tank dimensions, the modified ship 


parameters became: 


J = 1269.02171 + 11.4155 
= 1300.176 metric tons-meters-sec 
K = 683.280 - 3.121 


Ul 


680.159 metric tons-meter 


After four iterations of the design procedure, the 


following final values were obtained: 


W = 120 meter 
a = 1.982 meter 
2 
A = 1.982 meter 
Z 
a = 0.454 meter 
R = G.0 meter 
D = - 1.0 meter 
a = 1300.236 metric tons-meters-sec 
3 
B = 59.746 metric tons-meters-sec 
5 
K = 680.185 metric tons-meters 
5 
a = 124.276 metric tons-meters-sec 
an = 5.722 metric tons-meters-sec 
5 = 65.012 metric tons-meters 
Z 
J . = 266.152 metric tons-meters-sec 
3 


These values were then substituted into the transfer 
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function, the roots were located, and the frequency response 


was obtained. 


®(s) _ Ono 23 ES + 0.0460 s + 0.523 
W(s) ss + 0.164 s + 1.485 s + 0.08575 + 0.441 
(73) 
zeroes at s = -0.0230 + j 0.723 
poles at s = -0.0135 + j 0.642 
and s = -0.0684 + j 1.031 


The frequency response is graphically illustrated in 
Figure 2-2. The DSL/360 program which generated the plot is 
listed in Appendix 8. The negative decibel low frequency 
response demonstrates the ability of the passive tank systen 
to overcone the natural tendency of the ship to ride normal 
to the waveslope, making it ride upright with respect to 
earth reference instead. 

Compacing this response with the response of the 
unmodifiei ship, given in Figure 1-2, it was observed that 
the effect of adding passive anti-roll tanks was to place a 
pair of complex zeroes at the ship's own natural frequency 
to reduce the resonant peak, from 24 db at 0.7 rad/sec, to 
two lower peaks of 14 db at 0.6 rad/fsec and 10db at 1.0 
rad/sec, and a notch of -11 dab at 0.7 rad/sec. In sea state 
6, this translates to maximum ship roll angles of 8.5 
degrees at 0.6 rad/se2c, and 15.8 degrees at 1.0 rad/sec, a 
Considerable improvement over the 34.8 degree maximum roll 
angle of the unmodified ship. 
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Figure 2.2. Frequency Response of Ship with Passive Tanks 
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ieee SHEP STABILIZATION BY ACTIVATED ANTI-ROLL TANKS 


A. MODIFIED EQUATIONS OF MOTION. 


Activating an anti-roll tank system, mathematically, 
only involves adding a forcing function term to the tank 
equation of motion. Physically, however, it can be 
accomplished in basically two different ways. These methods 
are pneumatic and hydraulic. A pneumatic effector operates 
by alternately creating an overpressure or a partial vacuum 
condition in the airspace above the tank-fluid surface in 
the vyertisal legs of the U-tube, while an hydraulic effector 
generally consists of an arrangement of pumps in the 
horizontal cross duct. The effector type chosen for this 
analysis was a constant speed controllable pitch propellor 
pump, so chosen for its linear blade pitch angle vs. flow 


Saatacteristic | Blagoveshchensky, 1962 |. 
1. Nomenclature. 


In addition to the nomenclature of the passive tank 
System, the following are the parameters and variables of an 
activated tank systean: 
Dependent variables of the forcing side: 
@ = instantaneous pump-bladie pitch angle (commanded) 
Parameters of the forcing variable: 
B = velocity head-loss coefficient of the pump 
ie static head coefficient of the pump 


Torques produced on the forcing side: 





Ka- B#6 = Counter-rolling torque due to pump 
Q 
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eaei@diticatiion of Parameters. 

As indicated by the expression for counter-rolling 
torgue, the addition of a pump to the tank system has the 
effect of increasing the overall damping coefficient of the 
tank-fluid. All other parameters of the homogeneous side 
remain unchanged. The pump should be chosen such that the 


composite tank damping factor is 0.45 | Vachanaratana, e730) 


B = Bt + B (74) 

t p 

i. 0.45 (75) 

B = 0.45 a (76) 
9 


Seeeeotuations of Motion. 


The equations of motion for the activated anti-roll 
tank system come from equation (23), with Q, = Ke : Q, = Kyo 
; q, = @~, and die 0 [ Chadwick and Klotter, 1953]. 


Ship: (J @¢+Bo+tK oo) + (I @6+K 6) =K yp (77) 
S S = St is S 


Tank: (J O9+K%)+ (9 6+B6+K 6) =Ka (78) 
st t = t t D 


Following the procedure used for the passive ship-tank 


system, @(s) was found to be 


$(s) = we (S)),-- ae) (79) 
6 (s) = (K SJ [s° tts Jo .)s + (Kk JJ |v (Ss) (80) 
1 Ss ieee co 
a 
_ 2 
(5) = (K /3 (I 73)s + (K /3_)] A (s) (81) 
A = as defined by equation (40) 


B. CONTROLLER DESIGN AND FREQUENCY RESPONSE. 


The problem of ship roll control is a frequency domain 
problem: controlling against a (time domain) stochastic 
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disturbance with a generally known (frequency domain) 
Spectrum. Therefore the classical approach to control was 
pursued, designing on the basis of feedback compensation. 
It was assumed that roll angle, rate, and acceleration 
Signals would be available from shipboard sensors (gyros and 
accelerometers). The controller transfer function was 
considered to be of the form 
C(s) A(S)/ @ (Ss) 


il 


=Ks + KS + K (82) 
Z 1 0 


For purposes of plock diagram development and analysis, 
the following block transfer functions were defined, based 


upon equations (79), (80), and (81): 


SS(s) = ae / v(s) 
2 
=K /J + (B /d + (K /J 
/I_[s + B/I)s + (KI) |72 
2 2 2 
= w [s +2fwstow a (83) 
n 1} n 
35S) = ace / A(s) 
2 
= K/J Ke flys + (Kk /J )|7 4 
Dae Sc Ss t s 
2 2 
= kK /J ks + |74 84 
Ae [ xs _ (84) 
SST(s) = ST(s) / SS(s) 
2 2 2 2 
= ae | Xs rw Vis + 2008 + . (85) 
G(s) = KJ [s* + (B/J )S + (KK /J )] 
ss ee, ee St 


ee ce ttt 
aa 
K /3_ E + (BL/T)s + (K 5.) ] 


2 2 2 
= wW ls t2hwstw | (86) 
een n 


n 
are ls +2twste | 


n 
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H(s) SST (s) C(S) 


4 3 2 2 2 
1/K [XC_s #XC_s +(XC tw ys tw cs tac | 
S 2 1 O- Sn. 2 nq n 0 


s + 2f$ws +w (87) 
fon n 
Where 
c = K K 
Z BieZ 
c = K K 
1 p 1 
C = K K 
0 p 0 


The “open loop" transfer function, then, was 


F(s) = G(s) / [1 + G(s)H(s)| 
2 
= (K /3_)[s Pasey )s + (kK fd )] Pa 
A (open) 
- 2 
A(open) = A- (K_/3_)[s + (B/I)s + (R73) | + 


i /da)) XC ; ne ; (XC +(K /J )C ) ; 
+ S + Ss + + S 
s [ Z 1 0 toa t *2 


K K 
* Rae 4. vs a oa 


(89) 
The overall, or "closed loop" response became 
T(s) = P(s) / w (s) 
Sac) 7) 01 + F (5) 
K /J ) 5 + (B /J )S + (K SJ ) (90) 
= Ss S 
eee M7, 


A (ciose 
losed At (1/9 [ xc exc ag (KGe+ (K fic on 
close = S S 
= ( 5) 2 1 0 ce tC 


+(K /J )C S #(K SJ )C a4 
(K/3 YC Ss +(K /T)C | (91) 
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The Characteristic Eguation of the closed loop response, 
A(closed), was 


4 3 Z 
s +As t+ AS + AS #48 = 0 (92) 
2 1 0 
A = res) BJ+t JIB + J C 
3 s t s t s t st 1 
2 (93) 
J J-+3 .C JJ +J Jd 
st st 2 st st Ss. t 
A = J J K J+ J K + BBe- 2K J 
2 $s t Sat s t Sant st 
(94) 
JIJdg*tJ C JJ +J 
Ste St 2 s t st 
Uwe. + Ke 
Sie 20 t7.2 
+ 
JJ 
See 
A = J J BK* KB + KC 
1 s t Ss t s t 
(95) 
JI J+J C J Gg =J 
Ss Ss Ss st s t 
A = J J K K - K + K C 
J s t s t ie 
(2c) 
Jeu fro 8 26S JI J -J 
st st Z Sateeeees ¢ st 


The block diagram development is illustrated in Figure 
3-1. The first design approach attempted was the modern 
control theory state variable feedback method. The 
difficulty with this method was the inability to define a 
tine-domain objective or cost functional. 

The next design approach attempted was with the Computer 
Aided Design of Optimal Compensators (CADOC) progran 
[ Anastassakis, 1977]. This method consisted of selecting a 
Frequency-domain cost functional and a desired frequency 
response curve, and used the function minimization procedure 
of M.J. Box. There were two difficulties encountered in 
applying this program to the problem of ship anti-roll 
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See irse state of block diagram reduction. 


Figure 3.1. Block Diagram Development for Activated 
Ship/Tank System 
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Mey Final block diagram. 





(d) Equivalent block diagram for servo theory analysis. 


Figure 3.1. Block Diagram Development for Activated 
Ship/Tank System 
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control. The first problem was that a portion of the 
feedback path in the physical System was fixed by 
architecture, while the CADOC program required all 
coefficients in the feedback path to be adjustable. The 
second problem involved selecting an attainable frequency 
response for CADOC input. It was concluded that the CADOC 
program was best suited for optimizing a close design, 
rather than for initiating a design. 

Another design approach attempted was in the root locus 
parameter space. Two separate procedures exist here, and 
both were attempted. The first procedure was the standard 
root locis approach, fixing two of the three adjustable 
gains and varying the third. This method amounted to an 
exhaustive search technigue, attempting to locate the best 
combination of the three gains to define system roots with 
adequate jamping. 

The second root Locus parameter space procedure involved 
uSing a parameter plane computer program, PARAM A, which 
exists at the Naval Postgraduate School and is available 
from Dr. S.J. Thaler, but has not been documented in any 
literature to his knowledge. This program varies two gains 
which appear in linear combinations in the coefficients of 
the system characteristic equation, for a fixed value of the 
third gain. The output of this program is a two-dimensional 
parameter plane showing families of curves of constant 
damping ratio and constant natural frequency as functions of 
the two gains. 

The problems associated with both root locus parameter 
Space procedures were the exhaustive search nature of a 
three-variable system, and the exclusions of the system 
closed-loop zeroes from the analysis. This was noted to be 
a problem when the gains found by these procedures were used 
in a freyuency response study of the activated ship/tank 
system, and large resonant peaks were observed. 
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After attempting all of the above methods and abandoning 
them for the reasons stated, an exhaustive search of 
closed-loop frequency response curves was conducted, using 
the program listed in Appendix C. An acceptable design was 
accomplished when C = 110, C = -150, and C = -80. This 
corresponded to placing the compensator roots at s = -2.4388 
and s = +0.5638, which improved the overall system damping 
and provided the frequency response illustrated in Figure 
3-2. This response was characterized by two low resonant 
Means, One at 0.582 rad/fsec of 2.85 db (linear gain 1.39), 
and one at 1.38 rad/sec of 1.96 db (linear gain 1.25). The 
Minimum response between the two resonances occured at 0.861 
rad/sec, and was -4.75 db (linear gain 0.58). 

In sea state 6, the maximum waveslopes at the resonant 
frequencies are 1.44 degrees at 0.582 rad/sec, and 8.11 
degrees at 1.38 rad/sec, and the maximum waveslope at 0.861 
rad/sec is 3.15 degrees. Using this controller design, the 
maximum roll angles that the ship should experience in sea 
state 6 would be 2.0 degrees at 0.582 rad/sec, 1.83 degrees 
at 0.861 rad/sec, and 10.1 degrees at 1.38 rad/sec. The 
wave heizht power spectral density at those frequencies, 
according to tne spectra of Pierson and Moskowitz, are 205 
Et“ /Hz at 0.582 rad/sec (0.093 Hz), 80 ft ° /Az at 0.861 
rad/fsec (0.137 Hz), and 12 £t*/Hz atepleso brad/sec (0.2129 
Hz). Therefore it was concluded that even though the pure 
mumbers indicated a poor. roll control at the higher 
frequency resonant peak, the waves at that frequency lack 
sufficient power to cause the ship to attain the maximum 


theoretical roll angle. 
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BpeeveNe | RESPONSE, ACTIVATED TANKS 
meeUoew EOOP RESPONSE, MAGIDB) VS. LOG W 
: / \ /\ 


eee, 
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Figure 3-2. Frequency Response of Activated Ship/Tank 
System, Proportional Plus Rate Plus Acceleration 
Comeco! 
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Ven ween C lS Orevesoaonk DISPLACEMENT VARIATIONS 


The analysis of ship motions and roll control thus far 
had been conducted under the assumption that the 
displacement of the vessel was constant. Over the short 
term, the assumption may be true enough, but for a_ ship 
underway, consuming fuel and stores, the assumption cannot 
hold true beyond a few days at best. 

AS a ship's displacement decreases, the underwater 
sucface acea of the hull decreases also, in keeping with the 
Principle of Archimedes. This in turn leads toa lowering 
of the Canter of Buoyancy, which results in a lowering of 
the Metacenter. The effect on the metacentric height, 
however, 1S not so clear, as the exact location of the 
Center of Sravity in the ship depends more on the placement 
of the onboard weight than on the total displacement. [In 
general, fuel tanks are located low in a ship; thus 
consumption of fuel is removal of low weight, which shifts 
the Center of Gravity upward, decreasing the metacentric 
height. It should be noted here that ballasting empty fuel 
tanks with sea water has the effect of restoring the ship as 
meee as PISSible to the original condition of statical 
Stability which existed at the beginning of the voyage, and 
extends the time frame of applicability of the assumption of 
constant displacement. 

Varying the ship's displacement therefore varies the 
matural frequency of the ship, and effectively de-tunes the 
anti-roll tank system. The only physical adjustment that 
can be made to a passive tank installation underway is to 
vary the quantity of water in the U-tube. This has’ the 
effect of changing the limits of integration of the design 
eguations, and therefore changes the parameters of the tank. 
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Additional adjustments are available on the activated 
tank system. These are the gain of the pump, which isa 
function of the pump rotational speed, and the feedback 
gains from the sensors to the pump controls. Thus, ifa 
method exists for determining on-line the various’ ship 
parameters as they change, the controller for the activated 
tank systen can be made adaptive. Various schemes may be 
considered, from Kalman filter type parameter estimators to 
FFI devices, and this field is left as a topic for further 


possible cesearch. 
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Vee ONCLUSIONS 


mie stuudying ship roll control, i1t was necessary to 
become familiar with the nature of ocean waves. The 
waveslope spectrum of the ocean which was generated, when 
used in conjunction with the Pierson-Moskowitz wave power 
Spectra, gave much insight into the problem. It was learned 
that the se2 surface does not and can not assume the slopes 
that an observer on a ship often estimates. 

A typical ship on the ocean magnifies the waveslopes 
near the resonant frequency by a factor of about 15. By 
installing passive anti-roll tanks, this waveslope 
Magnification factor was split into two resonant peaks, one 
Mren a factor of 5, and the other with a factor of 3. 
Activating the tank system further decreased the waveslope 
magnification factors to 1.4 and 1.3, respectively. 

Several assumptions have been made over the years in the 
study of roll control by U-tube tank systems, and some of 
them were sources of problems in the design of physical 
systems. The first assumption was that a ship has constant 
parameters: displacement, metacentric height, and natural 
period. This assumption holds true only if the ship does not 
coosume any fuel or stores, or if these weights are replaced 
With appropriat2] quantities of suitably placed ballast. 

Another assumption was that the tanks could be placed in 
the ship in such fasShion as to drive the mutual-coupling 
Mmmecktia coefficient to zero. This assumption had _ two 
limitations. The first pitfall was the direct connection 
between this assumption and the assumption of constant ship 
parameters. The second pitfall was that the tanks had to be 
placed in an impossible location on the ship to accomplish 
the objective of the assumption. 


De 





Roll control by pump activated tanks is perhaps the best 
form of roll control for slow moving vessels, or for vessels 
conducting operations either at slow speed or stopped, where 
activated fin stabilizers lose their effectiveness. i 
snould be noted that a pump activated tank system can he 
used to cause rolling, as an aid for breaking out of ice. 
It can also be used to correct list, including damage-caused 
list, as a safer alternative to counter-flooding. Therefore 
the pump activated U-tube tank system is the most versatile 
form of coll control that a multi-purpose vessel can have 
installed. 
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APPENDIX A 


De eyepa@el i etal esponse of 


The Fcequency Response program was written in Digital 
Simulation Language (DSL/360) to be run on the IBM 
System/360 machine at the W. R. Church Computer Center, 
Naval Postgraduate School. To use the program, a _ deck of 
plot identification cards must be appended following the 
FORTRAN Subroutine, in accordance with the OSL/360 user's 
manual. 

The labels and variable names used in the program are 
defined in the Glossary for Frequency Response Programs, 


Appendix D of this thesis. 
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APPENDIX B 


2rogram Listing for Frequency poe of 
See with Passive Tanks, with 
ORTRAN Subroutine CALC 


The Frequency Response program was written in Digital 
Simulation Language (DSL/360) fon be run On the IBM 
System/360 machine at the W. R. Church Computer Center, 
Naval Postgraduate School. To use the program, a deck of 
plot identification cards must be? appended following the 
FORTRAN Subroutine, in accordance with the DSL/360 user's 
manual. 

The labels and variable names used in the program are 
defined in the Glossary for Frequency Response Programs, 


Appendix D of this thesis. 
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Meo grameoistindg £fOF Fre Boueue: Bea ous’ of 
Activated Ship/Tan em, With 
FORTRAN Subroutine CALCA 


The Frequency Response program was written in Digital 
Simulation Language (DSL/360) Ope run On the IBM 
System/360 machine at the W. R. Church Computer Center, 
Naval Postgraduate School. To use the program, a deck of 
plot identification cards must be appended following the 
FORTRAN Subroutine, in accordance with the DSL/360 user's 
manual. 

The labels and variable names used in the program are 
defined in the Glossary for Frequency Response Programs, 
Appendix D of this thesis. 
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APPENDIX D 


Glossary for Frequency Response Programs 
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Ao, A1,A2,A3 
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DL 

CLIOP 

eon P 


GIMAG 
GREAL 


IM 
IMA 
JS 
od 
JT 
KS 
Rod 
LOGW 
MAG 
MASH 
MDB 
MDBH 


OLOOP 
Pat 


MAIN 


MEANING 

Coefficients of characteristic equation 
Damping coefficient of ship 

Overall damping coefficient of tank-fluid 
Closed-loop transfer function 

Compensator (controller) transfer function 
Denominator polynomial of open-loop transfer 
buses Lon 

Open-loop transfer function 

Imaginary part of open-loop transfer function 
Real part of open-loop transfer function 
Closed-loop transfer function 

Imaginary part of open-loop transfer function 
Imaginary part of open-loop transfer function 
Roll inertia coefficient of ship 
Mutual-coupling inertia coefficient 
Pauveteont anertia coefficient of tank-fluid 
Seaecre Tiaguting cOGCrricient Of Ship 
Equivalent static coefficient of tank-fluid 
Natural logarithm of radian frequency 
Magnitude of open-loop transfer function 
Magnitude of closed-loop transfer function 
Decibel gain of open-loop transfer function 
Decibel gain of closed-loop transfer function 
Numerator polynomial of transfer function 
Open-loop transfer function 


Phase of open-loop transfer function 
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SYMBOL 
A 
AS 


BEAN 
Bie 
BS 
Br 
Bert 


DABS 
DCAL 


Dass P 


E1,E2 


GM 
GRAY 


ieee, 13,14 


JS 


Denominator polynomial of open-loop transfer 
PUnet Lon 
LaPlace transform variable 


Ship-tank response in feedback path 


SUBROUTINES 


MEANING 

Cross-sectional area of vertical leg of U-tube 
Cross-sectional area of horizontal leg of 
U=uibie 

Transverse dimension of ship (beam) 

Velocity head-loss coefficient of pump 

Damping coefficient of ship 

Overall damping coefficient of tank-fluid 
Sener Lercion coefficient of tank 
fluid 


Perpendicular distance from center of rotation 


damping 


of ship to horizontal cross-leg of U-tube 
Absolute value (magnitude) of DCAL 
Calculated 


Mmutual-coupling inertia coefficient to zero 


value Que D tou LOrce the 
Displacement of ship 


Incremental diffeneance between current and 


previously calculated values ym ship 
parameters 
Metacentric height 


Acceleration of gravity 


Height of fluid column in vertical leg of 
U-tube 
Evaluation of integrals as described in main 


body of thesis 


Roll inertia coefficient of ship 


oy 





JS1 
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KRAD 
KS 
Kot 
KT 


LEN 
LMDA 


RHO 
RHIG 
TJ 
suse 
TPHI 
WIDE 
WN2 


Roll inertia coefficient of unmodified ship 
Mutual-coupling inertia coefficiant 

Equivalent inertia coefficient of tank-fluid 
Radius of gyration of the mass of the ship 
Starlomnlgnting GOerEIcrent Of Ship 

Statice Lrilgneing coOetnicient Of unmodified ship 
Equivalent static coefficient of tank~-fluid 
Longitudinal dimension of a vertical leg of 
the U-tube 

Longitudinal dimension of ship (length) 

Design ratio A 

Number of iterations of design procedure 

Lever arm of U-tube 

Mass density of tank-fluid 

Weight density of tank-fluid 

Temporary Js 

Temporary KS 

Natural period of roll of ship 

Transverse dimension of vertical leg of U-tube 
Natural frequency of ship-tank system, squared 


Design ratio X 
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APPENDIX §& 


Waveslope Spectrum of the Ocean 


The waveslope spectrum of the ocean was generated from 
the relationships of the trochoidal wave model. tS 
intended fO be used an conjunction with the 
Pierson-Moskowitz wave power spectra. 

This appendix includes the listing of the DSL/360 
program which generated the spectrum, and listings of the 


spectra for sea states 3, 4, 5, and 6. 
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